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Abstract 

Short range meson-exchange mechanisms such as p — > ujtt contribute 
significantly to the amplitude for pp — > ppn near threshold in addition to 
pion rescattering. The uncertainty of the latter contribution, which depends 
on poorly determined parameters, is estimated. Allowing for this uncertainty 
and taking the short range meson exchange contributions into account makes 
it possible to reproduce the measured cross section. 
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The reaction pp — > ppir near threshold is peculiarly intriguing because 
the relatively large value of the momentum transfer (« 360 MeV/c) implies 
a very small impulse approximation term (Fig. la). Furthermore, the pion 
rescattering contribution (Fig. lb) is strongly suppressed because the dom- 
inant low energy pion-nucleon (ttN) scattering amplitude is the isovector 
Weinberg- Tomozawa term, which cannot contribute to this reaction |1| . The 
isoscalar pion rescattering term can be evaluated using chiral perturbation 
theory (ChPT), which relates this contribution to the off-shell 7rA scattering 
amplitude. Previous calculations ||| based on a 0(q 3 ) determination of 
this amplitude from the isospin-even scattering length and the sigma-term M 
found that the rescattering term leads to a significant destructive interference 
with the impulse approximation term. Part of this destructive interference 
arises from intermediate A(1232)'s (Fig. lc) ||. The net result is a near 
cancellation between the pion rescattering contribution based on the chiral 
Lagrangian and the amplitude for single-nucleon pion production (this can- 
cellation does not appear in off-shell extrapolations unconstrained by ChPT 
0). As a consequence, the empirical cross section |6|, [7|, ^ has to be un- 
derstood in terms of non-pionic two-body mechanisms of shorter range than 
that associated with pion rescattering. 



That shorter-ranged mechanisms should play a significant role in the re- 
action pp — > ppn° was already implied by the demonstration that the single- 
nucleon and the pion-exchange amplitude as conventionally calculated with 
the (now ruled out) on-shell approximation for the 7rA scattering amplitude 
predict a cross section smaller than the measured one by about a factor of 
5 [§, [nj. Subsequently it was found that the effective scalar and vector ex- 



change components of the nucleon-nucleon (AA) interaction imply relativis- 
tic two-nucleon operators of short range (Fig. Id) that contribute significantly 
to the cross section near threshold |TT| , |12|| . The same mechanisms, through 
PCAC, had earlier been shown to provide a natural explanation JT3] , |H| of 
the observed large enhancement in heavy nuclei of the axial charge of the nu- 
cleon |T5|, [To] . The cancellation between the single-nucleon and pion-exchange 



mechanisms shows, however, that the short-range operators implied by the 
scalar and vector components of the AA interaction cannot explain all of 
the measured cross section for pp — > ppir . 



We consider here meson exchange mechanisms in which the pion is pro- 
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duced at a transition vertex that involves two different heavy mesons (Fig. le): 
most relevant are pu, ttV and rjao exchanges. We also consider the effect of 
the intermediate iV(1440) resonance that can be excited by exchanged a and 
lo mesons (Fig. If). These contributions have a sizeable uncertainty range 
because of the values (and relative signs) of the various coupling constants 
are not well known. More dramatic, however, is the theoretical uncertainty 
associated with the pion rescattering contribution (Fig. lb), which is esti- 
mated by a re-evaluation of the rescattering contribution with another set of 
parameters: a 0(q 2 ) determination from ttN sub-threshold parameters [[l? 



In this case, the interference with the one-body term, while still destruc- 
tive, is far less complete. Finally, the uncertainty regarding the NN wave 
functions at small separations is estimated by using two different realistic, 



modern NN potentials, Argonne V18 [18] and Reid93 [19]. 



It turns out that the pu and, to a lesser extent, ir'a and r]a exchange 
mechanisms provide sufficiently large matrix elements for the reaction pp — > 
ppir so that, when combined with the short range exchange mechanisms 
that are associated with the NN interaction, they can provide an - at least 
qualitative - explanation of the empirical cross section. 



Consider first the pio exchange contribution to the two-nucleon pion pro- 
duction amplitude (Fig. le), as derived in Ref. pOfl . The key ingredient in 
this amplitude is the npu vertex, which was considered by Gell-Mann, Sharp 



and Wagner [21] in the p-dominance model for the decay to — > ixixix: 



r 



9-k 



e^sd^p'LU^n 



Here p v and lo\ are the Proca field operators for the p and uj mesons and 7? 
is the pion field. The irpuj vertex is related to the triangle anomaly, as illus- 
trated explicitly e.g. in Ref. ||22|| , and hence it can be viewed as fairly well 



established. The numerical value for the coupling constant g npuJ is ~ —10 
[Em (the phase is fixed as in Ref. p3|). 



The discussion of the poo exchange pion production amplitude is com- 
pleted by noting that we employ the usual uNN and pNN interactions with 
vector and tensor coupling. For the u- and p-nucleon coupling constants we 
use the values g 2 p /^Tx = 1.2, k p = 6.1 and g^/^Tc = 25, = respectively, 
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where, as usual, k p {w) represents the relative strength of the tensor coupling. 
These values correspond to the Bonn boson-exchange model for the NN in- 
teraction [23J and imply that the p and u exchanges are treated as effective 
meson exchanges rather than pure meson exchange interactions. We also 
assume that g p g^ > 0, in line with the quark model. The pNN and uNN 
vertex form factors are denoted F p (k) and F^ik) respectively. Their expres- 
sions may be inferred from boson exchange models for the NN interaction. 



The complete pu exchange amplitude for neutral pion production near 
threshold is evaluated to be 

T p,uj _ _.gpgug*pu ^ q F p (k)F w (k) 



4m 2 N m„ (k 2 + m 2 p )(k 2 + m 2 ) 
[(2 + k p + kJ(A? 2 S -P-k-PE-k) 
-2z(l + k p )(1 + K u )&ffM x <r (2) ■ k\. (2) 

Here we follow the notation of Ref. [0: uo 2 = q 2 + m 2 is the energy of the 

(on-shell) pion produced with momentum q in the center of mass; p (p') is 
the center of mass momentum of the incoming (outgoing) proton labelled "1" 
(those of proton "2" are opposite); k = p — p' is the momentum transferred; 
P = p + p'; <jW is the spin of proton i; and S = — a^. In Ref. |2(| only 



the numerically most important last term in this expression was retained. 



The puj exchange operator of Eq. (2) is insignificant for nuclear pion ab- 
sorption rates in near threshold, but it is of comparable magnitude to the 
(small) isospin symmetric pion exchange contribution to the pion production 
(or absorption) amplitude. As a consequence it gives a numerically signifi- 
cant contribution to the matrix element for the reaction pp — > ppn . 

In view of the dominant ttt] decay mode of the ao(980) meson and the 
important role the effective isospin 1 scalar exchange mechanism plays in the 
iViV interaction, the r/ao exchange mechanism (Fig. le) could be expected to 
be of similar significance for nuclear pion production near threshold as the pui 
exchange mechanism. The Trrjao vertex can be described by the Lagrangian 

£7^0 = mao9 * vao [mla ■ ttt] + r nr)ao a ■ d^d^r]}, (3) 
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where 3q is the isovector ao(980) meson field and rj the pseudoscalar ?7-meson 
field. Here g nV a is the non-derivative coupling constant and r nvao is the ratio 
between two-derivative and non-derivative couplings; higher-derivative inter- 
actions could be considered as well. We take the decay width T(a — > ni]) 
of the ao(980) to be 175 MeV, which is the mean of the empirical range 50- 
300 MeV pB| . We emphasize that the measured value of T fixes only one 
combination of the two parameters g wmo and r nr}ao . The derivative term is 
important in the decay width because of the large three-momenta of the pion 
and eta, but less important near threshold where the three-momentum of the 
pion vanishes. Hence the contribution of the wqao vertex to ir° production is 
very sensitive to the parameter r nr)ao . For example, with r^ao — 0, one finds 
I^tt^o! — 3.5 which would lead to a huge contribution. The values and 3.5 
are similar to those found in the linear basis of 't Hooft's extension of the 
linear a- mo del, which explains the U(l) anomaly EEL and would generate a 



large contribution to the pion production amplitude. There is, however, no 
justification for choosing = 0, as the rest of the calculation employs a 
pion field that couples via derivatives in the chiral limit. Rather, naturalness 
implies [rv^J ~ 1; e.g., in the appropriate non-linear basis, 't Hooft's model 
gives r nvao = —2, while assumed meson saturation of ChPT theory parame- 
ters in the meson sector gives [^7j r nr)ao ~ —0.75. Assuming [r^J ~ 1, the 
decay width is controlled by the two-derivative term, and |r wr)ao ^ w | ~ 0.21. 

The rjNN and a Q NN coupling Lagrangians are taken to have the usual 
pseudovector and scalar forms. The coupling constant f V NN is determined 



as fr)N n = { m 'r l /^ r nN)griNNi where g V NN = 5 J24|. Again interpreting the 
ao exchange mechanism as an effective representation of the isospin 1 scalar 
exchange component of the iViV interaction, we take the coupling constant 
g ao NN to have the value 5.9 indicated by the Bonn model |2l}] . The 7/iViV 
and a NN vertex form factors are denoted F v (k) and F ao (k) respectively, 
and can be obtained from the same reference. 

The complete rjao exchange contribution to the tt° production amplitude 
is then determined to be 



f F ri(k)F ao (k) F1 t ^ u, 

m 



2 



T V,a _ -i-^f 7}NN g aoNN g 7: - v a ° [1 + r^—JL-jE . k, (4) 

(m* + k 2 )(m* o + k 2 ) 



in the same notation as before. Obviously neither the sign nor the magnitude 
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of this contribution is determined by the decay width alone. In the follow- 
ing, this contribution is estimated by assuming the magnitude of the square 
bracket to be 1, as suggested by naturalness, and that the overall sign is such 
that this amplitude interferes constructively with the impulse approximation. 

One may proceed in a similar way to incorporate exchange of heavier 
mesons that can decay into 7r° and another neutral meson. There is, however, 
only scant information on the coupling of such heavier mesons to nucleons. 
An exception are the isovector mesons around 1.3 GeV: the ai(1260) and 
the 7r'(1300). Both have large decay widths, and could be expected to decay 
predominantly into 7r(7r7r)s_ wave . Arguments based on chiral symmetry and 
Regge behaviour can e.g. yield a significant decay width for ai(1260) — > 
ire, where e is a broad scalar more or less degenerate with the p ||28|| . An 
analysis that included both ai(1260) and 7r'(1300) has yielded a very 
small branching ratio for this process, however, while giving r(7r' — > Tie) = 
345 MeV. To be definite, we will concentrate on the 7r'(1300), and use this 
decay width as input. We consider the Lagrangian 

£w = ^j^lmlTr' ■ txo + r^d^il' ■ d^na], (5) 

where ff' is the isovector 7r'(1300) meson field and a the effective scalar field. 
Again, here g^-ov is the non-derivative coupling constant and r nan i is the ratio 
between two-derivative and non-derivative couplings; the pion decay constant 
f n was introduced for normalization. And again, assuming Ir^^l ~ 1, the 
decay width is controlled by the two-derivative term, and \r na . n ig WCT7r /\ ~ 0.7. 



The 7r'(1300) is not included in most boson-exchange NN potentials, but 
these have typically a large ttNN form factor parameter A„- ~ 1.3 — 1.5 GeV. 
However, there is a considerable body of evidence, see e.g. |30| and references 



therein, that this form factor is actually much softer. In particular, the use of 
hard pion-nucleon form factors, typical in nuclear physics, in calculations of 
deep inelastic structure functions leads to much larger breaking of the SU(3) 
flavor symmetry of the q sea than is observed |31], |32| . A value of 730 ± 100 
MeV is preferred by those analyses. 



The reduction of the effects of one-pion exchange in NN scattering caused 
by such a soft form factor may be compensated for by the effects of the ex- 
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change of a 7r'(1300) [Q. In particular, a form factor parameter A^ = 800 
MeV can be enforced in a one-boson-exchange potential without loss of fitting 
quality if the 7r'(1300) is introduced simultaneously, coupling to the nucleon 
as an ordinary pseudoscalar, with a coupling constant g\, NN jAm = 100 [[33 



In the same model, the a has its usual coupling to the nucleon, with a cou- 
pling constant g 2 aNN l^ = 8.38. The respective form factors |33j are denoted 
here F v (k) and F ao (k). 



In the same notation as before, the corresponding amplitude (see Fig. le) 



is 



, . ml F a (k)FAk) ^ q . r 

2m N f w (ml + k 2 ){ml, + k 2 ) 2m^ 

(6) 

Similarly to the rjao case, we assume in the following that the square bracket 
is of 0(1) and indeed is taken as unity in our calculations. We also assume 
that g a NNgnNN > as in the linear a- model, and that g^NNg^NN < and 
g-Ktjir'go-wK < as in quark models. 

As the AT(1440) is the lowest | + nucleon resonance it is expected to con- 
tribute to the S'-wave pion production amplitude in a similar way to the 
NN pair contributions that are implied by the Lorentz-invariant structure 
of the AW interaction. As the effective pion rescattering vertices include 
all ttN s-channel resonances, only amplitudes that arise from excitation of 
virtual A/"(1440) resonances by short range exchange mechanisms should be 
considered as separate two-nucleon pion production amplitudes. The most 
important of these are a- (or (^)s^ waye ) and cu-meson exchange (Fig. If). 

The appropriate vertices for derivation of the amplitude for the contri- 
bution that involves excitation of intermediate A/"(1440)'s through a- and 



cu-exchange are given in Ref. |34J. The aNN coupling constant g a is approx- 



imately given by the Bonn- model for the AW interaction [23]. The coupling 
constant g* for the excitation of the A^(1440) resonance is estimated from the 
partial width for the N(UA0) -> A^(7r7r)fr° ave in Ref. EM to be g* ~ 1.1. The 



value for the uNN* coupling constant, <?*, is highly uncertain. In Ref. |54j 
it was estimated to be g* ~ 1.7 under the assumption that g^/g u — g^/go- 
The 7rAW(1440) coupling is described by a pseudovector Lagrangian. The 
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ttNN* pseudovector coupling constant /* is obtained as /* = 0.62 from the 
empirical decay width for the iV(1440) — > Ntt decay p4]| . We assume that 
the form factors for the aN(14A0)N and AT (1440) AT vertices are the same 
as the aNN and uNN vertices, F a (k) and F w (k), respectively. 



The iV(1440) amplitude is then 

rpR-a,ui = , f: u r I Q.fMk) Sm N - m* g^F^jk) \ g p 
m N (m*-m N ) q \ k 2 + ml m* + m N k 2 + m 2 J 

+ A; 2 + m 2 , 

The form of this operator is similar to that of the effective meson-exchange 
term that involves an intermediate NN pair. The ratio of the iV(1440) 
scalar meson-exchange part of Eq. (7) to the scalar exchange pair contri- 
bution is (g*/ g a )(f*f n / 'gA^ir) (2mAr/(m* — m^)) ~ 0.25. This shows that 
the A/"(1440) contribution enhances the scalar meson-exchange contribution 
to the S'-wave pion production amplitude by ~ 25%. The two vector ex- 
change terms in Eq. (7) give contributions of different sign to the matrix 
element for pp — > ppn . Their ratios to the corresponding u pair terms are 
-(gt/gu>)(fnU/gAm n ) (2m N /(m*-m N ))({3m N -m*)/(m* + m N )) ~ -0.15 
for the first and (g*/ guj){f^fn/ gAm-w) {2m^/{m* + m^)) ~ 0.05 for the sec- 
ond. This partial cancellation between the u exchange terms that involve an 
intermediate A/" (1440) ensures that the u exchange contribution to the cross 
section for pp — > ppn near threshold will not be very different from the one 



obtained in Ref. Ill 



cor- 



All the above contributions, together with the pair terms of Ref. ]TT| 
respond to contact ttNNNN interactions in ChPT. According to the modi- 
fied chiral power counting developed in Ref. 0], the pair and A^(1440) terms 
contribute to the leading order of such contact terms, while the exchange 
mechanisms contribute to higher orders, i.e. they are down by one power 
of m % . In that reference, the cross section for the reaction pp — > ppir was 
calculated on the basis of a reaction amplitude formed primarily of (1) the 
single nucleon pion production operator, (2) the pion exchange contribution 
as constrained by chiral perturbation theory, and (3) the short range contri- 
butions associated with intermediate A^A" pair terms. It was shown that the 
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resulting cross section is too small by a factor 5 — 10 because of the remarkable 
cancellation between the single nucleon and pion rescattering contributions. 
This cancellation depends crucially on parameters determined by nN scat- 
tering. 

We shall here be mainly concerned with the modifications caused by (i) 
the additional short-range contributions discussed above and (ii) changes in 
the numerical values of the parameters that affect the pion rescattering con- 
tribution. 



The numerical values of the different contributions of the graphs of Fig. 1 
are shown in Fig. 2. These are computed using the Reid93 potential to 
provide the initial 3 Pq and final 1 5'o wave functions. The quantities J refer 
to the amplitudes after some common kinematic factors and complex phases 
are removed; see Ref. [[| for a precise definition. The label impulse refers 
to the contribution of Fig. la. There are two pion rescattering curves. They 
both include the same A excitation contribution (Fig. lc) and a contribution 
(Fig. lb) proportional to 



4Cl+ M 




c 



2 m N 



where q are input ttN parameters, and Sm^ ~ 3 MeV is the strong in- 
teraction contribution to the neutron-proton mass difference. The two pion 
rescattering curves differ in the q's and C. The one labelled cl (for chiral 
loops) uses the value of C = —2.31 obtained in Ref. M and corresponds to 
the 0(q 3 ) parameters q as given in Ref. 0. It is obvious that the contribu- 
tion of the impulse term is almost entirely cancelled by this pion rescattering 
term. This is why the corresponding cross section in Ref. |fj] was very close to 
zero without including the effects of the nucleon-pair terms associated with a 
and io exchange. The other rescattering curve (labelled st for sub-threshold 
expansion) corresponds to the 0(q 2 ) parameters q as given in Ref. |17 



which give a much smaller C = —0.29. (The subthreshold expansion makes 
use of an analytic description of the nN with numerical values determined 
by data |3"5j| .) As a consequence the cancellation is less complete in this case. 
Both of these curves were obtained with a soft pion form factor; employing a 
hard form factor increases the magnitude of these results by less than 10%. 
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This is because of the effects of NN repulsion in the final state 1 S'o pp wave 
function and the centrifugal barrier in the initial 3 Po wave function. The 
integrand for the pion rescattering contribution peaks at about 1.3 fm. 



The other curves in Fig. 2 display the shorter-range effects (Fig. ld,e,f). 
The effect of the new pun has the same sign as the impulse term and is 
comparable to the change in the rescattering term; it is also of similar size as 
the a pair term. The new r)ao7r° and an'it meson exchange terms are smaller 
under the (conservative) assumptions of this paper. The curve labelled u> + o 
includes the effects of the pair diagrams. The effects of the excitation of 
the iV(1440) are also shown. Including the iV(1440) enhances the sigma 
exchange term by about 25%, but the omega exchange term hardly at all 
because of the cancellation noted above. Therefore the effect of the iV(1440) 
is small. The results in Fig. 2, do have a wide uncertainty margin because 
of the uncertain meson couplings. In particular, using a value of 7v CT7r / much 
smaller than unity would increase the value of the related contribution by a 



big factor. Likewise, a larger value for g* |36[ would also increase the cross 
section. 



We also note that the magnitude of these shorter-than-pion-range terms 
is consistent with expectations from the modified chiral power counting of 
Ref. @. Although the convergence of ChPT here is not that impressive, 
there is no evidence of failure either. 



The cross sections, computed using all of the mechanisms of Fig. 1 as well 
as the very small recoil term of Ref. || are shown in Fig. 3, together with 
experimental data from Refs. |6|, [7], §]. The full curves represent the complete 
result for the Reid93 potential |19j, the lower curve using the 0(q 3 ) rescatter- 
ing parameters, the higher the 0(q 2 ) parameters. Also shown are the dashed 
curves with the corresponding results for the Argonne V18 potential JT^]. 
For a given potential, the two different set of rescattering parameters give an 
estimate of the uncertainty generated by the rescattering mechanism. For a 
given set of rescattering parameters, the two potentials give an estimate of 
the uncertainty generated by the nuclear short-range dynamics. The overall 
theoretical uncertainty from these sources is quite large; it arises because the 
cross section is basically very small. The inputs (both pp wave functions and 
various ttN scattering coefficients) are not known well enough to provide a 
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more accurate calculation. There is also additional uncertainty in the input 
coefficients of the meson-exchange Lagrangians of Eqs. (3) and (5). Further- 
more, two-pion exchange diagrams can be expected to contribute appreciably, 
but have so far not been calculated. Within this large inherent theoretical 
uncertainty, meson exchange can explain the observed cross section. 



Further development of ChPT should lead to a reduction of the uncer- 
tainties associated with the irN scattering coefficients, pp wave functions, 
and the mesonic coupling constants. A calculation to order (q 3 ) and com- 
parison with the sub-threshold expansion of Ref. |35j should lead to better 
description of 7rN scattering data. If the chiral potential of Ordonez et al 
[ |T7j ] were developed to yield a good description of NN data up to the ener- 
gies relevant here a consistent treatment of the amplitude for 7r° production 
and the nuclear wave functions based on the same general chiral Lagrangian 
could be achieved. Finally an improved treatment of the meson-exchange 
diagrams is possible by a chiral calculation of the relevant terms as higher 
order loop diagrams. 



The principal conclusion is that the influence of the pun - and to a 
lesser extent the other shorter-than-pion-range mechanisms considered here 
- is considerable. This idea can be used to understand independent experi- 
ments. In particular, the PCAC relation between the axial charge operator 



of the nuclear system and the S-wave soft-pion production amplitude flTl 
and the present results imply that the axial charge operator should have a 
significant pu axial exchange current term. As this will contribute to the 



nuclear enhancement of the effective axial charge of the nucleon [13|], this 
suggests a possibility for further reducing the remaining uncertainties in the 
magnitude of the short range mechanisms by data on first forbidden nuclear 



/^-transitions [113, 16] 
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Figure captions 

Fig. 1 Contributions to the amplitude for pp — > ppn . (a) Impulse 
approximation, (b) pion rescattering, (c) pion rescattering through an inter- 
mediate A, (d) pair graphs, (e) meson exchange contributions, (f) excitation 
of an intermediate iV(1440) resonance. 

Fig. 2 Individual contributions to the amplitude for pp — > ppir as 
a function of p', the relative momentum between the two final protons for 
laboratory kinetic energy of 300MeV. The Reid93 potential is used. 

Fig. 3 Computed cross sections for two different NN potential mod- 
els. All of the terms of Fig. 1 are included as well as that of the recoil term 
ofRef. [2]. 
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